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Abstract 

An expression for the internal energy of a fluid element in a magnetized, anisotropic plasma is 
derived from first principles. The result is a function of entropy, particle density and magnetic 
field, and as such plays the role of a thermodynamic potential: it determines in principle all 
thermodynamic properties of the fluid element. In particular it provides equations of state for the 
magnetized plasma. The derivation uses familiar fluid equations, a few elements of kinetic theory, 
the MHD version of Faraday's law, and certain familiar stability and regularity conditions. 



PACS numbers: 



I. INTRODUCTION 



A. Objective 

The equilibrium states of thermodynamics depend generally on a separation of time scales: 
they are equilibria only if processes involving some longer time scale are ignored. Thus one 
studies the vapor-pressure equilibrium of a glass of water without concern for the fact that 
the glass itself will eventually evaporate. 

Separation of time scales commonly pertains to the physics of magnetized plasmas, allow- 
ing thermodynamic ideas to provide useful illumination. However the full panoply of ther- 
modynamics, including the calculation of thermodynamic potentials, is rarely employed. In 
fact it is sometimes said that thermodynamics applies only to the fully equilibrated plasma, 
where the confining field becomes irrelevant: the eventual equilibrium state of a plasma 
immersed in a magnetic field is affected by that field only through intrinsic spin of the 
charged particles pQ. This stringent perspective misses the point that, at small collision fre- 
quency, there can be sufficient time-scale separation to speak usefully of magnetized plasma 
equilibria, in which the field plays an important role [2]. 

Here we consider the example of a magnetized, anisotropic plasma — a plasma in which the 
pressure tensor shows distinct parallel and perpendicular components [3]. We intentionally 
ignore the fact that, on some longer time scale, Coulomb collisions will erode the anisotropy. 
We consider a single fluid element in such a plasma, which is allowed to interact with 
neighboring elements and with the magnetic field B. These interactions allow the element 
to perform work on its environment, and the first task of our study will be to understand the 
work performed by an anisotropic fluid element. For simplicity we consider a single plasma 
species, and ultimately suppress the weak collisional interaction between species. 

Our final result is an expression for the internal energy of a magnetized, anisotropic 
plasma in terms of its natural variables: the entropy, the particle density and the magnetic 
field. Expressed in this way the internal energy becomes a true thermodynamic potential, 
from which other potentials, such as the Helmholtz free energy, are easily found. These 
potentials provide all the available thermodynamic information about the plasma system 
[I]. In particular the thermodynamic potentials readily provide equations of state. 

By confining attention to a single fluid element, we restrict the analysis to local thermo- 
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dynamics. Thus questions concerning the global equilibrium, such as the stability of fluid 
profiles, or the configuration of the confining magnetic field, are outside the purview of this 
work. We find however, that even within the local framework interesting and non-trivial 
thermodynamic conditions are revealed. 

An early application of thermodynamics to magnetized plasma is due to Fowler [5] (see 
also [6]), who used estimates of plasma free energy to derive approximate bounds on insta- 
bility growth rates and fluctuation levels. More recent thermodynamic analyses of plasma 
fluctuations [7J [S] have depended upon kinetic calculations of the perturbed entropy. Ther- 
modynamic calculations have proven similarly useful in the study of dusty plasmas; see, for 
example, the work of Avinash j9]. 



B. Organization 



The first four sections of this paper derive from fluid equations an expression for the re- 
versible work performed by a fluid element in a magnetized, anisotropic plasma. Elementary 
thermodynamics and kinetic theory provide, from the expression for the reversible work, a 
set of partial differential equations for the internal energy as a function of entropy, density, 
and magnetic field: Eqs. (30), (31), and (32). In Section [V] we find the general solution to 
these differential equations. After imposing appropriate regularity and stability conditions, 



we obtain the main result of this work, Eq. (58). This equation presents the internal energy 
as a thermodynamic potential in its natural variables. We explore some of its consequences, 
including the corresponding equations of state, which are expressed in terms of two different 
sets of variables. 



II. FLUID CONSERVATION LAWS 



Thermodynamic aspects of fluids have been treated in many works, notably [TOl [TTJ . Our 
presentation of this section differs from these classical works in two ways: we emphasize 
the notion of fluid element and we leave the pure thermodynamic setting by appealing to 
general constraints that arise by assuming an underlying single particle species plasma with 
no internal degrees of freedom. The results of this section set the stage for the results of 



Sec. Ill, where we further generalize by considering the role of the magnetic field. 
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A. Energy 



Within a purely thermodynamic/fluid mechanical setting one has the following general 
energy equation [TUHT2"] : 

3v 

— + V-(uV + q)+p:VV = W, (1) 

at 

where u is an internal energy density that is assumed to govern the thermodynamics of the 
fluid, p is a stress tensor, p = trace(p)/3 is a pressure, V is the Eulerian fluid velocity, q 
is a heat flux density, and W is the external power (rate of doing work), including energy 
exchange with other energy sources that might be present, which for our application could 
come about by interaction with another plasma species. We will leave W unspecified. 

Assuming a microscopic theory consisting of any single particle-species plasma the fol- 
lowing equation can be obtained for general such kinetic theories [32]: 

3 -^ + 3 -pV-V + p:VV + V-q = W (2) 



V ■ V 



where 

d _ d 
~dt'~ dt 

is the usual convective (material) derivative. 
Consistency of Eq. ([l]) with §2§ gives 

u= \p- ( 3 ) 

If we couple Eq. ^ with the assumption that our fluid is in local thermodynamic equilibrium 
as described by the energy representation, and assume u is only a function the entropy and 
volume, then pressure is given by differentiation with respect to volume (see e.g. [H ITT] ) 

P = n-, (4) 

where ti := u/n and n is the particle density. Equations ^ and Q immediately imply 

u = An 5/3 , (5) 

where A depends only on entropy; i.e., we obtain the thermodynamic internal energy function 
for an adiabatic monatomic gas. 
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Our development of Sec. Ill where we treat anisotropic magnetized plasma, is a general- 
ization of this basic idea, where the thermodynamics is generalized to include the anisotropic 
effect of the magnetic field (cf. Eqs. (30), (31), and (32) below that are analogous to Q 



and (|3j)). In the remainder of this section we develop several notions that will elucidate our 
approach and be useful for later analysis. 

The term p : VV in the above equations represents work done by the fluid stress. It is 
helpful to express this work in terms of two traceless tensors: the rate of strain tensor 



and the viscosity tensor 



The result is 



U aP = i {d a V p + d p V a ) - l -5 a f7 ■ V 



7T a /3 = Pa/S -$ a pP- 



nil ) 

— + -uV ■V + ir:U + V-q = W. (6) 

at 3 



B. Entropy 



We denote the entropy density by s(x, t), the entropy flux density by s(x, t) and the local 
rate of entropy production by Q(x,t). Thus we have 



Thermodynamics prescribes the flux 

.s = sV 4- 

T 



8 = sV+±, (7) 



where T is the temperature. It follows that 

£ + .v.v + Iv., = e + £.H. ( s) 

The use of the thermodynamic relation ^ becomes questionable if the distribution function 
does not resemble, in some approximation, a moving Maxwellian. Thus at this point we 
implicitly assume that resemblance. A more detailed discussion of the distribution and of 
^ is presented in subsection III A 



After solving ([8]) for V • q and substituting the result into (j6]), we find that 

du ^ds /5 

T h - 

dt dt \ 3 



/ " T ?+ (-u-Ts) V-V + 7z:U + q-V\ogT + Te = W. (9) 



The leading terms in this equation have a simple thermodynamic interpretation, which we 
consider next. 



C. First law for a fluid element 

We consider the physical system consisting of a fluid element, small on the scale of 
plasma gradients but containing many particles. The environment for the system is the 
surrounding plasma. The element is defined by the particles it contains, and moves with 
those particles; therefore its population N is fixed and the chemical potential will not appear 
in our development. However, the fluid element volume 

V = N/n 

changes according to 

dV = -Vd logra. 

Here n is the plasma density. Since the energy of the fluid element is U — uV we see that 
an energy change dU is given by 

dll = d(uV) = V(du — udlogn). 

Similarly the entropy change dS of the fluid element is given by 

dS = V(ds — sd\ogn) . 

We use these formulae to compute 

dU -TdS = Vdt 



du ds dlogn 

— T- (u — Ts) — : — 

dt dt v 1 dt 



The first law of thermodynamics states that the left-hand side of this relation is the reversible 
work 8W r performed on the elemental system, 

dU - TdS = 5W r = Vdt w r , (10) 

where w r is the rate of change of work-density. Thus the thermodynamic law is expressed 



as 



du ds , ^ . d log n . . 
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After noting that 

^ = -v-v, ( 12 ) 



we substitute (11) into d9j) and find that 



TQ + w r = W-q- VlogT- \uV ■ V - ir : U. (13) 

Evidently every term on the right-hand side of this relation must describe either reversible 
work (w r ) or irreversible dissipation (TO). In some cases the categorization is obvious; for 
example, 

-\uV-V = p ^P^- = -{Vdty l pdV (14) 
reproduces the reversible work done on an ideal fluid. Similarly collisional heat-conduction, 

q c = -KVlogT, 

contributes to — q • V log T a positive-definite term 

-q c - VlogT = (^VlogT) - VlogT 

that obviously belongs to TO. But for other terms the identification is not obvious. The 
most interesting term is that involving the plasma viscosity. 



III. MAGNETIZED PLASMA VISCOSITY 
A. Small-gyroradius decomposition 

Notice that the magnetic field has not entered the formalism explicitly up to this point. 
It does so through the form of the plasma viscosity, which we now consider. The viscosity is 
computed using a small gyro-radius ordering, so at this point we depart from general theory 
and specialize to a magnetized plasma. 

In typical contexts the particle distribution function f(x,v,t) for a magnetized plasma 
has the form [HI M, El HE] 

/ = fu + /a + f g (15) 

Here the first term denotes a Maxwellian distribution, centered at the mean flow velocity V; 
the second term is a correction to the Maxwellian, independent of gyro-phase, that includes 
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the stress anisotropy Ap = p\\ — p±; and the third term, which depends on gyro-phase, 
describes gyration about the magnetic field. Both correction terms are first-order in the 
small gyro-radius parameter 

5 = pr/L 

where px is the thermal gyro-radius and L is a typical scale length for system gradients. 
The distribution will in general contain second- and higher-order terms, but they have no 
effect on the present analysis. 



Equation (15) has two well-known (see, for example |14j ) consequences. First, it confirms 



the thermodynamic form of the entropy flow, given by ([7]); straightforward calculation from 



(15) shows that this form remains valid through first order in S. The second consequence 
concerns the form of the generalized viscosity 7r; one finds that the viscosity of a magnetized 
plasma decomposes into three parts: 

Here the first, gyrotropic term is that emphasized by Chew, Goldberger and Low (CGL) [3]: 

1 



izgt ee Ap \^bb - -I 

with b = B/B and Ap = p» — p±; the second term is conventional collisional viscosity, 
discussed in the following subsection; and 7r gv represents gyroviscosity. The detailed form 
of the gyroviscosity tensor is well-known [T31 E], and has been used in many works (e.g. 
[TBI HE]) but it is not needed here. For in fact 

7z gv : U = . (16) 

In other words gyroviscosity, while having important effects on momentum evolution, does 
not enter plasma thermodynamics. 



Equation (16) pertains to the gyro- viscosity tensor as it is given in most the literature; 



see, for example, [131 Q31 [T7] . Other contributions to gyroviscosity [T5] may contradict (16); 
such contributions would yield an additional term, not considered here, to the reversible 
work. 
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B. Collisional viscosity 

We comment on the form of the collisional viscosity since conventional treatments [TJ1 
[T?j [T8"] are coordinate-dependent. Any second-rank tensor W can be decomposed into 
components aligned with the magnetic field, 

W\\ = b(b ■ W) + (b ■ W)b - bb(b -W-b), 

and perpendicular components Wj_ = W — W\\. We use this decomposition to write the 
standard collisional viscosity formulae as [HI [TFJ [19] as 

TTc = 7T C || + 7T C _|_ , 

with 

Tc|| = —fJ'oUw, 7T c j_ = —fJ,lU± . 

Here the two viscosity coefficients /j,q and fix differ. We note that a distinct parsing, in which 
part of the collisional viscosity is lumped with what we call ir gt , is also possible [H]. Our 
iTgt evidently refers to stress anisotropy alone. 

We note that this decomposition is also useful with regard to gyroviscosity [20J. 



C. Viscous work 



We combine the above formulae to compute the viscous work (strictly, viscous power 
density) specific to a magnetized, anisotropic plasma: 



7T : U = noU\\ : U\\ + fJnU± : U± - Apb ■ U ■ b . 



(17) 



The first two terms here are positive definite and obviously contribute only to 0. But the 
last term can have either sign; this is the gyrotropic work, which we denote by 



w 



at 



-Ap b-U b = -Ap 



b-(V||V)-^V-V 



(18) 



For a physical understanding of gyrotropic work, we consider the case of uniform magnetic 
field and incompressible flow, in which 



gt 



-ApV||Vf|. 
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The factor ViiVji corresponds to contraction (or expansion) along the direction of the mag- 
netic field, for which the relevant force is p». But to preserve the volume, this distortion 
must be accompanied by an opposite change in the directions transverse to the field. Since 
this second change acts against the force p±, and since it enters with opposite sign, the work 
done must be proportional to Ap = p\\ —p±. It is clear that this work, like the ideal version 
—pdV, can be reversible. 



In summary, (13) has become 
TQ + w r 



W -q-V\ogT + Q v -pV -V 



Ap 



b-(V||V)--V-V 



(19) 



where 



e„ :=/i C/|| : U\\+mU x : C7± 
gives the collisional viscous dissipation. We associate the first three terms on the right-hand 



side of (19) with entropy production, and the remaining terms with reversible work: 

-pV ■ V - Ap 



UK. 



An alternative version is 



1. 



w,. 



b-(V||V)--V-V 



- P± V ■ V- Apb- (ViiV). 



(20) 



(21) 



IV. MHD VERSION 



A. Field evolution 



We need only one characteristic of MHD, the relation 

Vx£ = -Vx(VxB) 



(22) 



This requirement is much weaker than the statement B 2 V = E x B; it allows in particular 
for diamagnetic flow, provided Vnp = 0. Combined with the parallel component of Faraday's 



law, (22) yields 



dlogB 
dt 



b-(ViiV)-V-V 



(23) 



We note parenthetically that (j 1 2[) and (|23|) together require 

b-(Vi|V). 



d\og(B/n) 



dt 



10 



B. First law for MHD 



The MHD expression for gyrotropic work follows from (18) and (23): 



v ,'dhgB 2„ Tr 



Finally, the full MHD work is found from (20) or (21): 



w. r 



-P||V -V - Ap 
d log n 



dlogB 



-P\\ 



dt 

dlogV 



Ap 



dt 
dlogB 



dt 



- Ap 



dt 

dlogB 
dt 



We now return to (10), which becomes 

5W r = Vdtw r 

or, since V = N/n, 



N'HWr 



\dV-Ap^-dB 
B 



P -\dn-^dB 
n z nB 



Finally, from the first law 



dU 



we infer the relations 



TdS + 5W r 
TdS + N 



dU 

dn 
dU_ 

dB 



n 



\dn - ^-dB 
nB 



N 



P\\ 



n- 



-N 



Ap 
nB 



(24) 



(25) 



(26) 
(27) 



(28) 
(29) 



as given by Morrison [21 J, who obtained them by requiring energy conservation in CGL 
theory for general thermodynamics, i.e., arbitrary U(n, s, B). 

It is convenient to express these relations in terms of the normalized energy U = U /N — an 
equivalent potential since N is fixed. We have 

dU p\\ 



dn n 2 ' 
dU _ _Ap 

dB " ~n~B 



(30) 
(31) 
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In addition to the thermodynamic expressions (31) we have the relation 

'3 



U 



n 



:P\\ - 



(32) 



which was given in the original CGL paper [3], and emerged from the single species 



anisotropic kinetic considerations there. Note, Eq. (32 ) can be inferred directly by comparing 



with the usual energy density expression u — p/i^y — l), where 7 = (d + 2)/d with d the num- 
ber of degrees of freedom. With anisotropy one would expect u = P\\/ — 1) + P±/(j± — 1), 
and upon choosing d = 1 for 7y and d = 2 for 7j_, one arrives directly at (|32]). We note, 
however, it would be wrong to assume p\\ t ± ~ n 3 ' 1 , as we shall see in Sec.|V| 



The thermodynamic relations (31) together with the kinetic result (|32h constrain the 



form of the internal energy function - this we turn to next. 



V. INTERNAL ENERGY EXPRESSIONS 



A. General internal energy 



Upon inserting (30) and (31) into (32) we obtain 



U = -nU n + BU B 



(33) 



a linear first order partial differential equation that has the following general solution ob- 
tained by integrating the characteristic equations [22J: 

U = n 2 / 3 f(S,B/n 2 / 3 ) . (34) 



where / is an arbitrary function. Equation (34) implies the following expression for the 
internal energy density: 

u = n 5 / 3 / (S, B/n 2 ' 3 ) . (35) 



From (34) we obtain the following pressure relations 

P\\ = |n 5/8 / - \Bnf 
P± = 2 -n^f+ l -Bnf. 



(36) 
(37) 



We note, that we cannot further specify U without adding more physics. However, we 



note that / is not entirely free; in Sec. VP we discuss physical constraints on it. It is worth 



noting that in the limit B — > one simply obtains form (34) the adiabatic monatomic gas 
result. 
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B. Entropy dependence 



Now we turn to constraints on the entropy dependence of our general internal energy 



function of (34). We imagine the situation where our fluid element has equilibrated to a 



single temperature, even though the magnetic field can sustain anisotropic pressure. We 
define temperature by appealing to kinetic theory, where it is defined, like all quantities of 
moment equations, in terms of moments of the particle distribution function. Temperature 
as usual measures mean kinetic energy, and in terms of our fluid moments, it takes the form 

nn V 2m 2, 



n 



3n 



-U 



(38) 



where p = (pn + 2p±)/3 is the total pressure. Equation (38) essentially expresses classical 
equipartition. On the other hand, thermodynamically temperature is determined by the 
usual relation in the energy representation 



T = U S . 



After equating (38) and (39), then inserting (34) we obtain 



U = n 2/Z e 2S/Z g ( B/n 2/3) 



(39) 



(40) 



where g is an arbitrary function. 



C. CGL internal energy 

In the original CGL paper [3], the following equations were given: 



U^)=0 and £(af^)=0. (41) 
dt \nBJ dt\ n 3 J ( J 



Because entropy is advected, these can be used to further restrict U, since expressions (41) 
will be true if 

^ = c ± (S) and Pj f = c {l (S) (42) 



for arbitrary functions c± and c\i. Upon returning to (37), assuming / = h(S)g(B/n 2 / 3 ), 
and setting x := B/n 2 ^ 3 , we obtain the equations 

d± = T = sh 9ix) + \ 9 ' {x) (43) 
d \\ = f = ~y9{x) - —g\x) (44) 
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with d± and d\\ arbitrary constants. The solution of Eqs. (44) is g = d±x + d\\/(2x 2 ) and, 



therefore, consistent with (41), the internal energy function is 



!/3 /i(<S) U 



5 



n 2/3 2 fi 2 
,2 



dj. B + 



2 5 2 

whence we obtain the following expressions for the pressures: 

,3 



(45) 
(46) 



P± 



n 
IP 

d±h(S) nB 



(47) 
(48) 



Alternatively, instead of (41) the following are sometimes proposed: 

d (v\\B 2 ' 



d 



, ( P ~4)=0 and ; 
dt \ n 5 J dt 



rv 







(49) 



Performing the analogous calculation for these expressions, will produce the same internal 



energy function of (|46|). The authors of [23] use these relations and obtain 

and 



oc n 



p± oc n 



if 



where 



7f 



1 + 



HB/Bq) 



ln(n/no) ' ^" ln(n/no) ' 

with uq and Bq being constant reference values. Note these satisfy 7jj ff + 2^ e ± 



eff 



3-2 



ln(B/B ) 



(50) 



(51) 



5. The 



pressure relations of (50) are equivalent to our (47) and (48). This way of writing them 



demonstrates that CGL theory does not have constant polytropic indices for the two pres- 
sures, but ones that can be interprerted as having spatial dependence through n and B; 
evidently, it would in general be wrong to assume p»± ~ n 3,1 . Although our thermodynamic 
formalism has a single temperature, one can define 



T± := p±/n cx B T\\ := p\\ oc (n/B)' 



(52) 



which might aid intuition. 



We conclude this subsection by noting that the procedure of Sec. |VB| can be used to 
select the function h. 
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D. General constraints: nonnegativity, extensivity, and stability 



The function / of (34) is not entirely arbitrary and is subject to usual constraints of 
thermodynamics. To address these, we rewrite out internal energy function in terms of 
standard thermodynamic variables appropriate to the energy representation; i.e., entropy 
S = SN, volume V = N/n, and internal energy U = NU, where N is the fixed total number 



of particles. Thus (34) becomes 

U (N, S, V,B) = N {N/V) 2/3 f (S/N, B (V/iV) 2/3 ) . (53) 

The first comment to make is that any suitable internal energy function should be non- 
negative, and this is an elementary requirement on the function / for the relevant ranges of 
its thermodynamic independent variables. 

Next, it is well-known that the energy representation is the natural extensive one: the 
extensive internal energy is written in terms of the extensive particle number (here constant), 



volume, and entropy. The extensive property is obvious from (53) since U is an Euler 
homogeneous function of degree one in these variables: U(cN, cS, cV, B) = cU (N, S, V, B). 
Note that B has not participated in this scaling. This is because B is an intensive variable 
and we have opted to use it rather than the total magnetic moment, which is the conventional 
extensive variable for magnetic systems. 

Lastly we require thermodynamic stability. In the energy representation, equilibrium lies 
at minimum energy. Convexity of U assures us that unphysical behavior, such as having the 
pressure drop upon compression, will be ruled out. If U were to only depend on V and S, 
then the following local stability conditions would be necessary: 

U vv > and U ss > (54) 

as well as the Hessian condition 

UvvUss ~ (Uvs) 2 > (55) 

These inequalities place constraints on the function / as well as similar conditions involving 
the B dependence. We will consider a particular case below. 
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E. Polynomial form 



For purposes of concreteness we now expand the function g of (40) as a polynomial in 
x = B/n 2 / 3 . A second-order form is appropriate, since the higher-order terms introduce 
unphysical singularities; thus we obtain 

(B B 2 \ 

a O + Cll ^73 + a2 ^i73 ) • ( 56 ) 

where the coefficients a« depend on entropy alone. 

Next we express the internal energy of the fluid element in terms of dimensionless measure 
of density and magnetic field. To this end we consider the entire plasma macro-system, 
of which the fluid element is a part. We suppose that this system is characterized by a 
minimum density value, n m , and a maximum value of the magnetic field, Bm] the convenient 
dimensionless variables are then 

h:=n/n m , B:=B/B M 
and we find it convenient to introduce 

X(n,B) := B/h 2 ' 3 . (57) 
With these variables and some scaling we obtain 

u(S, n, B) = c (S)h 5 ^ [1 + ci(5)A + c 2 (S)X 2 ] . 
Here Cq evidently has the dimensions of energy density, while C\ and c 2 are dimensionless. 



Including the entropy dependence as described in Sec. |VB| and selecting the overall constant 
by taking the unmagnetized limit, gives the following: 

u(S, n, B) = u e 2S/3 n 5/3 (l + ai \ - a 2 X 2 ) , (58) 

where the constants ai j2 are fixed in each fluid element and 



with M being the particle mass. The expression of (58) can be viewed akin to the virial 



expansion for correction of the ideal gas law, correction here due to anisotropy. 
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Presently we consider constraints on the constants a\ t2 as well as the reason for minus 



sign in front of a 2 , but before doing so we observe that the pressures corresponding to (58) 



are given by (30) and (31): 



P|| = \u^ s l^ (1 + a 2 \ 2 ) 
Ap = -M e 2lS/3 n 5/3 5(a!A - 2a 2 X 2 ) . 



(59) 
(60) 



Now we apply the constraints discussed in Sec. VP to the weak field internal energy of 



(58). First we show that one can always choose the coefficients a« to guarantee that the 



energy density is non-negative for every fluid element in the plasma system. Our normal- 



izations guarantee that the quantity A of (57) satisfies 



A < 1 



(61) 



for every fluid element. It follows that the energy density will be non- negative provided 
A c > 1, where A c is the parameter value at which u vanishes: 



1 + a\\ c — a 2 \ 



0. 



It can be seen that A c = 1 only at a 2 = 1 + «i so the positivity requirement is simply 



ai > , < a 2 < 1 + a 



i • 



(62) 



We observe from (60) that this constraint allows both positive and negative values of the 



pressure anisotropy Ap. 

Next we consider stability. A system is thermodynamically stable if its thermodynamic 
potentials are convex functions of the intensive variables [3]. Thus we obtain the local 
stability criteria 

U m <0, U BB <0 



Straightforward calculation shows that 



nJAnn = --u e 2S / 3 (l + 5a 2 h- 4 / 3 B 2 ) 
y 



U 



BB 



2u e 2S / 3 a 2 n-V 3 



(63) 
(64) 



Hence thermodynamic stability introduces no additional constraints. Of course the ther- 
modynamic stability of a single fluid element does not by any means insure overall plasma 
stability. 
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F. Helmholtz free energy 



The Helmholtz free energy J-(T, n, B) is related to U by a Legendre transformation: 



J= = U-TS 



where S has been expressed in terms of the natural variables (T,n,B). Straightforward 



manipulation of (40) yields 



2 



Similarly, for the special case of (58) manipulation yields 

'2u ^n,B)h 2 ^ 



1 + lOE 



in terms of the abbreviation 



(65) 



(66) 



£(n, B) = 1 + ctiA — a^A 2 ■ 



The B = version of (66) is well-known. 



We differentiate to compute the entropy, 

S -- 



^ = ^ l0g (l 



and the pressures 



Th T 

p\\ = n 2 7 n = —(1 + a 2 \ 2 ) 
3nT 

Ap = -nBT B = - — (aiA - 2a 2 A 2 ) 



(67) 
(68) 



These equations of state restate (59) and (60) in terms of the variables (T, n, B) 



VI. SUMMARY 



We have found the functional form of the internal energy U of a fluid element in a 
magnetized, anisotropic plasma, as a function of its natural variables, the entropy S, the 



elemental volume V and the local magnetic field B. This function, given by (53 ), constitutes 



a thermodynamic potential, containing all thermodynamic information about the element. 
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In order to derive it, we first used standard fluid equations and the MHD version of Faraday's 
law to find how the fluid element performs reversible work. We then combined that result 
with standard thermodynamic laws to derive partial differential equations determining the 
dependence of U on the density n (which takes the place of V) and B. The general solution 
to those equations revealed that the magnetic field enters thermodynamics only through the 
combination B/n 2 ^ 3 - a quantity having the dimensions of magnetic flux. 

We next used the equipartition theorem to determine the entropy dependence of the 
internal energy. The resulting form generalizes the well-known expression for the thermody- 
namic potential of an ideal gas, and leads to new equations of state in which the magnetic 
field and anisotropy make key contributions. Thus we have found, in particular, how the 
adiabatic law used in conventional MHD is modified by anisotropy. 



We also a display a special case of the general result, Eq. (58), having quadratic form. In 
this case the thermodynamic potential involves two undetermined constant parameters. A 
simple constraint on these parameters simultaneously guarantees positivity of the internal 
energy and thermodynamic stability of the fluid element. 

The thermodynamic potential has experimentally testable consequences, including equa- 
tions of state. In particular it predicts a change in sign of the anisotropy m — p± as the 
quantity B/n 2 ^ 3 increases. Laboratory tests of such predictions would require measurements 
of the plasma on time-scales shorter than the time for collisions to relax the anisotropy [2] 



It would also be interesting to see if an equation of state built from (34) could provide a 
good fit to measurements of naturally occurring plasmas (e.g. [25]). 
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